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Abstract. We have evaluated the spatial distribution of 
dust in disks supplied by colliding or evaporating bodies 
spread in a small belt. The gradient of the distribution 
of dust expelled by radiation pressure is generally steeper 
than the one observed around /3Pictoris. It follows that, to 
generate the /? Pictoris disk distribution in such a way, one 
should extract relatively more small particles in the pro- 
duction process: thanks to radiation forces, these smaller 
particles have large eccentricities and thus could be seen at 
very large distances from their injection place. The evapo- 
ration process of comet-like bodies moving slowly inwards 
may yield the necessary particle size distribution. 

It is interesting to note that this dynamical description 
is able to account for the main observational properties of 
the (3 Pictoris disk, such as the power law distribution 
with possibly a change in slope, the asymmetry at large 
distances, and the total mass, all of which have remained 
unexplained up to now. If confirmed, this scenario may 
indicate that the f3 Pictoris disk could be looked on as a 
gigantic multi-cometary tail with its natural constituents: 
gas and dust. 
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1. Introduction 

Since the discovery of the dusty disk around (3 Pictoris 
(Aumann 1984, Smith & Terrile 1984), it has been the sub- 
ject of extensive investigations with all possible observa- 
tional techniques: visible and IR imagery (Smith & Terrile 
1987, Paresce & Burrows 1987, Lecavelier des Etangs et 
al. 1993, Gohmowski et al. 1993, Lagage & Pantin 1994), 
spectroscopy (Vidal-Madjar et al. 1994, Vidal-Madjar & 
Ferlet 1995 and references therein) and recently photo- 
metric variations (Lecavelier des Etangs et al. 1995). In 
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particular, the stellar light scattered by dust particles has 
been well observed by several authors for many years to 
follow a power law distribution with a slope in the range 
-3.6 (Artymowicz et al. 1989, Lecavelier des Etangs et al. 
1993) -4.3 (Smith & Terrile 1984). However, a convinc- 
ing explanation for this scattered light distribution is still 
needed. 

The dust distribution in the f3 Pictoris disk may be 
a crucial issue for understanding its origin. We present 
here a precise dynamical model which will be compared 
with the observations for the first time in a quantitative 
manner. 

Another key point revealed by the spectroscopic stud- 
ies of the gaseous counterpart of the /3 Pictoris disk is the 
existence of orbiting kilometer size bodies which some- 
times fall onto the star, giving rise to the observed spec- 
troscopic redshifted signatures due to ejected gas, and un- 
doubtedly to dust particles. Already in 1984, Weissman 
(1984) questioned whether the material around Vega-like 
stars is of cometary or asteroidal origin and suggested that 
particles could be continually supplied by sublimation or 
from collisions between larger bodies. 

The time scale evaluation over which dust particles 
are eliminated by the Poynting-Robertson effect or colli- 
sions (Backman & Paresce 1993) leads to durations much 
shorter than the estimated age of /^Pictoris 2 • 10^ years 
according to Paresce, 1991). This simple comparison sug- 
gests convincingly that the /SPictoris disk is not a remnant 
of planetary formation but on the contrary must be contin- 
ually replenished by secondary sources like evaporation or 
collision of small bodies. This seems in fact a very natural 
hypothesis (see for example Zuckerman & Becklin 1993): 
the presence of small bodies is a consequence of planetary 
systems formation (Lissauer 1993). 

Furthermore, CO absorptions have been recently ob- 
served towards (3 Pictoris (Vidal-Madjar et al. 1994) with 
the Hubble Space Telescope. The very presence of CO may 
also require a permanent source provided by evaporation 
of comet-like bodies. 
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Here we present a new argument: This replenishment is 

able to reproduce the main characteristics of the /3Pictoris 
dust disk, and in particular the spatial distribution of dust 
at large distances. In detail, it gives possible explanations 
for the following unexplained issues: 

1. The gradient of the scattered light follows a rela- 
tively well-known but unexplained power law (e.g. 
Golimowski et al. 1993). 

2. The distribution at large distances is obviously not 
axisymmetric (Smith & Tcrrile 1987, Kallas & Jewitt 
1995). 

3. The central part of the disk is relatively clear of dust 
(Backman et al. 1992, Lagage & Pantin 1994). 

4. The disk seems to be a "wedge" disk: the thickness in- 
creases with radius (Backman & Paresce 1993, Lecave- 
lier des Etangs et al. 1993). 

5. The slope of the scattered light distribution changes 
abruptly at about 100 AU from the star. If confirmed, 
this fact remains unexplained (Artymowicz et al. 1990, 
Golimowski et al. 1993). 

Furthermore, the possible explanation of these obser- 
vational facts could give new understanding to the follow- 
ing questions: 

— a) What is the mass of the dust disk, since when ex- 
trapolated towards infinite distances the disk mass di- 
verges (Artymowicz et al. 1989)? 

— b) Are there connections between the dust and the gas 
disks ? 

— c) Can the asymmetry in the observed dust disk be 
connected with the asymmetry in the longitude of peri- 
astron of the comets in the Falling-Evaporating-Bodics 
(FEB) model proposed to explain the redshifted spec- 
troscopic events (e.g. Beust et al. 1991)? 

Therefore, we will consider here a model in which the 

disk is replenished by a group of kilometer-size bodies 
(Section 2). Numerical calculations of the spatial distri- 
bution of dust are given in Section 3. We shall discuss in 
Section 4 the /? Pictoris disk in such a scenario and sum- 
marize the results in Section 5. 

2. The model 

2.1. Basic concept 

We assume that the dust is produced by cometary or 
asteroid-like bodies which create particles through mutual 
collisions or evaporation processes. The main assumption 
is thus that these particles have small initial velocities 
relative to the parent bodies: less than 1 km/s like in 
cometary production (Gombosi et al. 1985, Sekanina 1987) 
or in "Chiron burst" (Luu & Jewitt 1990). Thus, these ve- 
locities can be considered as negligible in comparison to 
the orbital velocities. In our model, this relative velocity 
is fixed at zero. However, as soon as a dust particle is 
ejected, it is perturbed by the radiation pressure; its orbit 



is different from the parent body one and tangent at the 

point where it was injected (Burns et al. 1979). For exam- 
ple, if a parent body in a circular orbit with a semi-major 
axis ao produces a particle with a ratio (3 of the radiation 
force to the gravitational force, the particle orbit has a 
semi-major axis 0^3 = oq (1 — /3)/(l — 2/3) and an eccentric- 
ity €0 = /3/(l — (3). The periastron is 0^3(1 — 6/3) = oq and 
the apoastron 0^3(1 + ep) = ao/(l — 2/?). 

Therefore, particles can be observed at distances from 
the central star much larger than the apoastron of the 
parent body. Thus, local perturbations on the distribution 
of the parent bodies like asymmetries, could produce ob- 
servable signatures on the dust distribution at very large 
distances. 

2.2. Particle size distribution 

Furthermore, because the parent bodies should not pro- 
duce single sized particles, we introduce in the calcida- 
tion a size distribution. For particles with radius s > l/z, 
p is correlated with the size of the particle by (3 ~ s~^ 
(Artymowicz 1988). We can assume a power law size dis- 
tribution as in Solar System cometary dust (Lien 1990 
and references therein), in collisionally replenished dust 
(Grecnberg & Nolan 1989), or in interplanetary medium 
dust (Le Sergeant & Lamy 1980). Thus, if we have a size 
distribution dn oc s'^ds, then we have dn oc fi^dfi with 
K = -q-2. 

2.3. Analytical consideration 
2.3.1. One point production 

First, let us consider a parent body in a circular orbit 
which generates a set of particles at a given point. Then 
the particles with (3 > 0.5 have hyperbolic orbits and are 
ejected from the system towards the interstellar medium. 
The other particles with /? < 0.5 follow orbits within a 
parabola which is the orbit of the particles with (3 = 0.5. 

One can evaluate the surface density of particles in the 
asymptotic direction of the parabola: for a given particle 
the true anomaly 6 follows the distribution law 

PeiO)de = -4^ — T^de (1) 

^ ' 27r(l-hecos6»)2 

In the asymptotic direction 6 = n, Pein) oc (1 — 

e)"^/2 cx r^/^{l + ao/ry/^. The probability that the par- 
ticle apoastron is between r and r + dr in the same di- 
rection is Prir) = l3^{dP/dr)dr. Since (3 = {1 - oo/r)/2, 
Pr{r) cx (1 — ao/r)^r~'^dr. Thus, taking into account the 
scattering cross section a oc /?~^ we obtain the surface 
density normal to the plane of the disk cr • rig: 

cr • 27rr • ns{r) oc r'^^^il - ao/r)^(l + ao/r)^/^ (2) 
Thus a ■ ns(r)r_>oo oc r~^'^. 
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2.3.2. Parent bodies in circular orbits. 

However, dust production can take place at every point in 

the orbit of the parent body. Here, we restrict ourselves 
to a parent body in circular orbit, thus, producing an ax- 
isymmetrical disk. Each particle has a probability //3(r) to 
be at a distance r: fpir) oc ((1 — ao/r)(2/3— 1 + 
The surface density normal to the plane of the disk is then: 



a ■ ns{r) oc 



1/2 



(r- 



(3) 



This surface density follows an r ^ law and is plotted in 
Fig. 1. 

Moreover, the inclination of the particles is the same as 

the inclination of the parent bodies. Therefore, the vertical 
distribution of the particles depends on the distribution of 
the inclination of the parent bodies. In the present model 
the thickness of the disk is increasing with radius and the 
volume density (n„) is proportional to the surface density 
[us) divided by the distance, i.e. n^(r) a r~^. 

It has to be noted that this conclusion is valid at dis- 
tances larger than the distance of the farthest parent bod- 
ies. Within this zone the opening angle must look smaller. 

From the ny{r) oc r""* law and Nakano's (1990) conclu- 
sion concerning the connection between the volume den- 
sity and the scattering hght distribution (F(r) oc n„ (r ) /r) , 
one can conclude that a belt of parent bodies in circular 
orbits should produce a "wedge" disk with a scattered 
light distribution F{r) oc r~^. 




2 4 
ln(r/aj 

Fig. 1. Plot of the surface density distribution calculated from 
Eq. 3. The abrupt density decrease at r ~ oo is simply due to 
the assumed single distance between the central star and all 
the parent bodies. It is thus an artificial consequence of an 
oversimplifiod hypothesis. 

For comparison, the hatched zone represents the distribu- 
tion range of slopes deduced by different observers around 
/3 Pictoris. 



3. Numerical results 

3.1. The Monte- Carlo simulation 

The analytical model of the previous section can only be 
solved for a limited number of parameters and parent body 
distributions. We have chosen to evaluate the dust dis- 
tribution by a Monte-Carlo method. We shall consider a 
given family of parent bodies, a dust production law and 
a particle size distribution, and we shall integrate the dis- 
tribution in the following way: 

For each of the N random particles, we first randomly 
select a parent body in the family and its position in its 
orbit; then we randomly apply the particle size distribu- 
tion and we calculate the position ri of the particle i after 
a time t where t is randomly chosen between and tmax- 
Finally, we evaluate F{r) the surface brightness along the 
midplane of the edge-on disk projected on the sky at the 
projected distance r. If Fi{r) is the contribution of the 
particle i, from Eq. 2 of Nakano (1990), we have 

Fi{r)ocr-^ r ni{A)ai{e)de (4) 
Jo 

where A is the distance to the central star and ai{9) 
is the scattering phase function. Nakano (1990) has al- 
ready shown that, except if the dust is much more steeply 



forward-scattering than in the Solar System, the isotropic 
scattering assumption does not change the observed slope 
of the scattered light distribution, thus ai{9) = fi^"^ with 
sin^ = r/A. The volume density is given by ni{h) = 
S{A — ri) ■ Pi/ri where Pi is the production rate of dust on 
the parent bodies, then Fi{r) oc J 5{A — ri)piH(A — 
r)rdA/(rif}fA^ A"^ — r^) where H is the Heavyside func- 
tion. 

Finally, the surface brightness at a given projected dis- 
tance r from the star is estimated by 

F(^ \-F(^ f nin -r).p^ 

Pir)=}_^F,{r)o^^^—= . (5) 

3.2. Results 

Results of runs with different initial conditions on the 
family of parent bodies are summarized in Tables 1 and 
2. The production rates per annulus of fixed width are 
given for each run. Several configurations have been tried. 
The first one corresponds to the analytical model of 
Sect. 2.3.2, we find again a = 5.0, which validates the 
Monte-Carlo model. In runs #2 and #3, we have tried dif- 
ferent distances and evolutionary times which show that 
the Poynting- Robertson effect does not change the results: 
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the accuracy for a is about 0.1. Prom #4 to #11, the par- 
ent bodies are in eccentric orbits with eccentricities uni- 
formly distributed between and 0.5 or between and 
1; it seems that bodies with eccentricities larger than 0.5 
do not contribute much to the distribution since a large 
fraction of ejected material is then on hyperbolic orbits. In 
runs #4 to #7, we have taken parent bodies with a pecu- 
liar longitude of periastron (a>), namely 155°. This value 
has been selected according to the Beust et al. (1991) mod- 
elisation of the Falling- Evaporating-Bodies, the so-called 
FEB scenario needed to reproduce the spectroscopic red- 
shifted signatures. In the four last runs of Table 1, uj was 
uniformly distributed between 0° and 360°. 

From all the results in Table 1, it can be seen that the 
slope a is always greater than or equal to 5.0. In Table 2, 
we test other parent body distributions in which the dust 
observed at large distances still is produced close to the 
star. Incidentally, another solution would be to assume 
that the disk is produced by a distribution of asteroids up 
to 1000 AU from the star. However this solution does not 
allow us to explain the observed asymmetries. 

Runs # 12 to # 15 take into account the Epstein gas 
drag with a gas density pgas — po(100 AU/r) and tem- 
perature Tgas = 20 K. Runs # 16 and # 17 are with dust 
size distribution characterized by = 5 — 10. 

4. The j3 Pictoris disk 

4..I. Orbiting- Evaporating- Bodies (OEB). Towards a solu- 
tion? 

The observed distribution in the external part of the 
/3 Pictoris disk {a < 4.3) is less steep than all calculated 
slopes in Tables 1. It has to be noted that a solution com- 
patible with the observed slope could be obtained with a 
gas density po ~ 10^ cm""^ at 100 AU. However, this is 
excluded by the upper limit deduced from spectroscopic 
or HI observations (Vidal-Madjar et al. 1986, Preudling et 
al. 1995). 

Runs #16 and #17 are also compatible with the 
dust distribution in Pictoris (Fig. 2). But they require 
dust size distributions {K = 5 — 10) very different from 
the expected ones if they are produced through colli- 
sions between parent bodies (Fujiwara 1979 and references 
therein). However they show us that we need a relatively 
larger amount of smaller particles which have larger apoas- 
trons. 

In the evaporation process of large bodies (like Chi- 
ron) , there is an upper limit s„iax to the size of the ejected 
particles extracted by gas (Luu & Jewitt 1990). Following 
Cowan & A'Hearn (1982) we have calculated this upper 
limit and the evaporation rate Z in molecules per second 
per unit area for CO and CO2 as a function of the dis- 
tance to /3 Pictoris, assuming L/jpic = SL© and a typical 
parent body radius Rbody = 100 km (Figs. 3 and 4). This 



upper limit is inversely proportional to the parent body 

size Rbody 

We have introduced this effect in our model. The re- 
sults are summarized in Table 3 for CO2; they would 
be similar for CO if the parent bodies were at larger 
distances. Here, we have taken the production rate 
p (X Z ■ R^g^y and a parent bodies size distribution 
dn = R'body^^'-^body wlicrc 7 = 3.5 (uncertainties de- 
duced from observations in the Solar System with 7 be- 
tween 3.2 (Whipple 1975, Hughes & Daniels 1982) and 3.8 
(Fernandez 1982) do not change the present results). 

In runs # 101 to # 103, we take a wide belt of parent 
bodies from 15 to 30 AU, without correlation between 
their radius Rbody and their distance from the central star. 
We obtain a ~ 4.8 — 4.9; the effect of cutoff in particle 
size is not efficient enough to explain the observed slope 
a -3.6 -4.3. 

However, if the evaporating bodies are coming from 
larger distances and their orbital parameters diffuse from 
a Kuiper belt-like zone towards a planetary-like zone, the 
smaller parent bodies lost their volatile material very early 
at larger distances whereas the larger bodies can evaporate 
downward at smaller distances. 

This scenario can be modelized in the following way: 
a parent body with a characteristic distance (which can 
be the periastron or the semi-major axis) r(t) = r^ — rt 
is considered as evaporating at r = r(T) only if the total 
mass of the previously evaporated gas is smaller than the 
available mass of volatile material. 

^ pZ{r{i))A^^Rl^^y\dt<^■p^Rl^^y (6) 

where A is the percentage of active surface of the body, 
^ the relative mass of CO2 in the parent body, p its density 
and n the CO2 molecular weight. Eq. 6 can be reduced to 

Rbody > Rref C ^(O^^V iZf ^(^^^ ^^^^^ Rr ef is the 

reference size of the smallest body which can evaporate 
at r = rref - For CO2 evaporation and rref — 20 AU, one 
obtain 

In runs # 104 to # 115, only those with Rref ~ 
30 — 40 km give the observed slope range. Moreover, due 
to the particle size cutoff, (3 > 0.4 there is an abrupt 
break in the slope at rbreak ~ «o/(l — '^0min) ~ 100 AU 
(see Fig. 5 for run #106). For CO, one can obtain r = 
6 • 10~^A/(5p)AU year~^, and there would be also an 
abrupt break in the slope since CO begins to evaporate 
from about 100-150 AU. 

4-2. Discussion 

We saw in the previous section that the observed gradient 
in the /? Pictoris disk can be explained by the present 
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Table 1. Slope a of the gradient of brightness in disks with different conditions for the parents bodies. In all these runs we 
assumed N=10000 and K=1.5. Poynting-Robertson drag is taken into account but there is no gas drag. For runs # 4 to # 7, 
the two slopes are for the both sides of the asymmetrical disk seen from the Earth. 



Run 




Parent Bodies Paxameters 






tmax 


Result: a 


# 


semi-major axis (AU) 


Eccentricity 




Production rate 


(years) 


F{r) oc r"" 


1 


20 


0. 


- 


p = const 


10000 


5.0 


2 


2 


0. 


- 


p = const 


10000 


4.9 


3 


20 


0. 


- 


p = const 


100000 


5.0 


4 


20-30 


0.-0.5 


155° 


p = const 


10000 


5.2-5.6 


5 


20-30 


O.-l. 


155° 


p = const 


10000 


5.2-5.7 


6 


20-30 


0.-0.5 


155° 


p (X r^^ 


10000 


5.3-5.5 


7 


20-30 


O.-l. 


155° 


p DC 


10000 


5.3-5.6 


8 


20-30 


0.-0.5 


0-360° 


p = const 


10000 


5.2 


9 


20-30 


O.-l. 


0-360° 


p = const 


10000 


5.3 


10 


20-30 


0.-0.5 


0-360° 


p oc 


10000 


5.2 


11 


20-30 


O.-l. 


0-360° 


p oc 


10000 


5.4 



Table 2. Same as the previous table, with different K; gas drag is taken into account. The parent bodies have semi-major axis 
o =20 AU and eccentricity e = 0. tmax=WOOO. For the runs # 16 the slope gradually changes from a = 3.7 at 20 AU to a = 4.7 
at 800 AU; and for run # 17 from 2.9 to 4.4 respectively. 



Run 


P distribution:K 


gas density at 100 AU 


Result: a 


# 


dn = 


(cm--') 


F{r) oc 7--" 


12 


1.5 


100. 


5.0 


13 


1.5 


1000. 


5.1 


14 


1.5 


10000. 


4.6 


15 


1.5 


100000. 


3.8 


16 


5. 


0. 


3.7-4.7 


17 


10. 


0. 


2.9-4.4 



model only if there is a larger amount of small particles. 

Now, there is an upper limit in the size of particles ejected 
by slowly evaporating bodies due to a balance between gas 
drag and gravitation. In fact, we obtain dust disks very 
similar to the /3 Pictoris one by assuming evaporation of 
kilometer size bodies which are diffusing from a Kuiper 
belt-like zone inward an evaporation zone. 

The time scale r ^ 10~^ AU year~^ in the diffusion 
of comet orbital elements, which is necessarily to explain 
P Pictoris's characteristics, is consistent with the few AU 
in 10 Myr deduced by Torbett & Smoluchowski (1990) 
and Levison's (1991) results on the Kuiper belt dynami- 
cal behavior between 30 and 40 AU around the Sun. The 
analysis of the evolution of the orbital elements of the 
parent bodies is outside the scope of this paper; more- 
over, the motion of such Kuiper belt-like bodies under the 
influence of a planetary system is probably chaotic (SchoU 
1979, Torbett & Smoluchowski 1990) and very dependent 
on this planetary system. The present model merely shows 
us that the dust distribution may be an important clue 
and may be explained by the present model if there is a 
long time scale perturbation of parent bodies by a plan- 
etary system which is compatible with the one predicted 
for the Solar System. 



Another evolution process for heating the bodies con- 
sists in the increase of the star luminosity (L) during 
its evolution along the main sequence. For a star of 
1.6 Mq, we have dL/Ldt = 2 • 10"^° year-\ if f is the 
distance of equal luminosity, this corresponds to f = 
f/2L • dL/dt. This solution would have the advantage of 
being "planetary system" independent. However, if the 
primary volatile is CO, f = 10~* AU year"'^; even if the 
volatile is CO2, f = 2.5 • 10~^ AU year~^. The time scales 
do not seem to be adequate. This mechanism would be ef- 
ficient only if we could have bodies with very large quanti- 
ties of volatiles and very small fractions of active surface. 

The assumption that CO or CO2 can be primary 

volatiles is not a problem since around an AV star they 
can evaporate at very large distances where the presence 
of young objects without a lag deposit of nonvolatile form- 
ing a crust is likely. On the contrary in the Solar System, 
we observe bright comets only when H2O can evaporate, 
simply because CO or CO2 evaporate only inside plane- 
tary regions where, for dynamical reasons, the presence of 
young objects is unlikely: Chiron seems to be an excep- 
tion but provides a useful analogy with what may happen 
around (3 Pictoris. 
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Table 3. Same as the previous tables, for particles ejected by evaporating gas. Here we assume that the gas is CO2 and 
p oc Z ■ Rlody Rref is the radius of the largest dead comets at rre/=20 AU. For runs # 105 to # 115 there is a slope break at 
r = rbreak, and both inward and outward slopes Uin and aout are given. For runs # 113 to # 115 the two given slopes Uin are 
for each side of the disk. We always took tmax=WOOO years. 



Run 




Parent Bodies Parameters 




Rfef 


"^break 


Result: 

CXin 


F{r) oc r-" 

OLout 


Brightness 
ratio 


# 


Periastron (AU) 


Eccentricity 




(km) 


(AU) 








101 


20-30 


0. 


- 


- 


- 


- 


4.8 


- 


102 


20-30 


0.-0.5 


0-360° 


- 


- 


- 


4.9 


- 


103 


20-30 


0.-0.5 


155° 


- 


- 


- 


4.9 


- 


104 


15-30 


0. 


- 


10. 


- 


- 


4.6 


- 


105 


15-30 


0. 


- 


20. 


50 


2.9 


4.5 


- 


106 


15-30 


0. 


- 


30. 


100 


2.7 


4.3 


- 


107 


15-30 


0. 




40. 


180 


2.5 


3.9 




108 


15-30 


0.-0.5 


0-360° 


10. 






4.6 




109 


15-30 


0.-0.5 


0-360° 


20. 


50 


2.8 


4.5 




110 


15-30 


0.-0.5 


0-360° 


30. 


130 


2.8 


4.3 




111 


15-30 


0.-0.5 


0-360° 


40. 


180 


2.6 


3.7 




112 


15-30 


0.-0.5 


155° 


10. 






4.6 


2.9 


113 


15-30 


0.-0.5 


155° 


20. 


50 


2.2-4.0 


4.5 


2.8 


114 


15-30 


0.-0.5 


155° 


30. 


110 


2.7-3.7 


4.3 


3.5 


115 


15-30 


0.-0.5 


155° 


40. 


180 


2.7-3.4 


3.7 


4.0 




Log J, r (AU) 

Fig. 2. Plot of the surface brightness (arbitrary unit) of the 
equatorial plane of an edge-on disk seen from the Earth cal- 
culated in run # 16. The slope varies from -3.7 to -4.7. For 
comparison, the dashed lines represent the power laws r~^''^ 
and r"* '^. 



As we shall see below, this model is also able to give 
answers to other important issues defined in Section 1. 



4-3. Asymmetry 

The distribution observed at large distances in (3 Pictoris 

disk is obviously not axisymnictric (Smith & Terrile 1987, 
Kallas & Jewitt 1995). However, relative keplerian mo- 
tions should remove such asymmetries: for two parti- 
cles on circular orbits with semi-major axis oi and 02, 
the time needed to put them in opposite side is t ~ 
0.4(aj ^^'^ — a2^^'^)~^ years if ai and 02 are in AU. With 
ai = 100 AU and a2 > 400 AU we have t < 450 years, an 
extremely short time compared to the age of the system. 
Thus, any asymmetry should quickly disappear. 

As noted above, it is however possible to find the 
asymmetry observed between the two extensions of the 
(3 Pictoris disk by assuming a peculiar parent body dis- 
tribution (see Fig. 6). Rims # 112 to # 115 have been 
obtained with a peculiar longitude of periastron of the 
parent bodies {lo = 155°). We found a brightness ratio be- 
tween the two extensions of the disk ranging from 2.8 to 
4.0. This result shows that the observed ratio of SW to NE 
extension brightness from 1.1 (Lecavelier des Etangs ct al. 
1993) to 1.5 (Kallas & Jewitt 1995) can be reproduced if 
there is a small additional proportion of bodies with pe- 
culiar longitudes of periastron like in the FEB scenario. 

However, there is another way to obtain the observed 
asymmetry, namely to assume that bodies with different 
longitudes of periastron move inwards with different ve- 
locities. For example, if bodies with a given longitude of 
periastron move more slowly inward, it produces a fainter 
disk extension with a steeper brightness gradient in one 
particular direction (Fig. 7). 
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Fig. 3. Plot of the evaporation rate Z as a function of the dis- 
tance from /JPictoris. The thick lines are for CO2 and the thin 

ones for CO. The evaporation is calculated for a steady-state 
energy balance and depends on the albedo Av of the parent 
body. The short-dashed lines are for A„=0.1, the long-dashed 
are for ^^=0.8 and the solid line are for Ay=0.5. In all the 
simulations we take ^,,=0.5. 
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Fig. 5. Plot of the surface brightness F{r) calculated from run 
# 106. There is an abrupt change in the slope at r = 100 AU. 
The slope is a = —2.7 inside and a = —4.3 outside this limit. 
For comparison, the dashed lines represent the power laws r~^'^ 
and r"* 
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Fig. 4. Plot of the radius of the largest particle which can 
be ejected from the surface of a parent body with radius 
Rbody=WO km and density Pbody = 1 as a function of /SPictoris 
distance. The thick linos arc for CO2 and the thin ones for CO. 
The different albedos are the same as in Fig. 3. 




Logio r (AU) 

Fig. 6. Plot of the surface brightness of the two extension 
calculated from run # 113. In this run all the parent bodies 
have a longitude of periastron of 155° relative to the line of 
sight. The brightness ratio between the two extensions is equal 
to 2.8. For comparison, the dashed lines represent the power 
laws r^^'^ and r^*'®. 
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Fig. 7. Plot of the surface brightness calculated from the addi- 
tion of runs # 108 and #115. In this combined run the parent 
bodies with a longitude of periastron of 155° go more slowly 
inwards. The brightness ratio between the two extensions is 
equal to 1.3. The fainter extension has a steeper slope and is 
represented by a long-dashed line. 

As in Fig. 1., the hatched zone represents the distribu- 
tion range of slopes deduced by different observers around 
B Pictoris. 



4.4- Mass of the disk 

The total mass of the /? Pictoris disk calculated from op- 
tical observations is 



M = 



ns(r, sj—pszTTrdrans 
o 



(8) 



where s is the particle size and dUg the size distri- 
bution. As discussed above we can take drig = s~^'^ds. 
Previously, authors assumed that s is independent of r 
the distance from the star. This assumption is the simpler 
one and others must be justified. 

This means however that: 



M oc 



J ns{r)rdr j s~^-^ds 



(9) 



The second integral diverges towards the largest bod- 
ies; the first one also diverges for rout — » oo if the slope 
of the scattering light is > —4 (that is the case according 
to Artymowicz et al. (1989), Lecavelier des Etangs et al. 
(1993) and for Gohmowski et al.'s (1993) NE extension). 
Thus to evaluate the mass, an arbitrary limit for the disk 
extension and particle size must be fixed. 

In the framework of our model both problems are di- 
rectly solved. At a given distance r, we have 2e < s < 



where ao is the distance of the parent body to the star and 
e is defined by e = sj3. We saw in Section 2.2 that e can 
be assumed to be constant. We obtain 



M oc 



M oc 



Aiip 



ns{r)2TTrdr 



(2er)/(r-Oo) 



-0.5 



ds 



167r- 



! roi 

-pV2l / 



(10) 



(11) 



This integral now converges and the computed mass is 
of the order of few lunar masses, i.e. of the same order as 
the lower limit evaluated by Artymowicz et al. (1989) who 
assumed an arbitrary outer limit of 500 AU and single-size 
particles of Ip. This mass is consistent with the total mass 
deduced from submillimeter observations (Zuckerman & 
Becklin 1993). 

4-5. Gas- dust ratio and CO detection. 

One can investigate the connection between the dust ob- 
served aroimd f3 Pictoris since 1984 and the gas disk with 
its stable component thanks to the recent detection of 
CO in UV lines with the HST (A^^^ « lO^^ cm'^, Vidal- 
Madjar et al. 1994). The very presence of CO may also 
need a permanent replenishment naturally provided by 
evaporation of comet-like bodies. 

We can evaluate the total mass of dust M4 which is 
associated with the evaporation of CO: 



Md = Mcoifi'dust/gasi'gas/CO 



td 

tco 



(12) 



2er 



where Mco is the total mass of CO in the disk, 
and ipgas/co are the mass ratios of dust to gas and gas 
to CO, (p is the mass ratio of the dust produced effec- 
tively kept in the disk and td and tco are the life time of 
dust and CO: due to photodissociation by UV interstel- 
lar radiation tco ^ 300 years (Vidal-Madjar ct al. 1994), 
and at a distance of 100 AU where CO begins to be va- 
porized, td ^ 10^ years (Backman & Paresce 1993). From 
the upper limit to particle size (/3 > 0.4), one can eval- 
uate (f ~ 0.1. With a simple geometry of a disk with 
an opening angle of 10", we can deduce from the ob- 
served CO column density a total CO mass of about 
Mco ~ 2 • 10^° kg. From observations in the Solar System, 
we can take ipdust/gas > 0.1 (Newburn & Spinrad 1989). 
Finally, i^gas/co depends on the composition of the evap- 
orated gas; it is between 1. if only CO is present and ~ 10. 
if all volatiles evaporate (Mumma et al. 1993). Thus, with 
'fii'dust/gas'^l-'gas/co ~ 0.1, Md is about onc hmar mass, 
that is of the same order as the total mass of the dust 
disk independently evaluated! 

This rough calculation shows that CO and dust are 
compatible with a common origin and produced by the 
same process: evaporation of comet-like bodies at large 
distances from the star. 
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5. Conclusion 

The model we have proposed is able to account for the 
main characteristics of the (3 Pictoris disk as natural con- 
sequences of the production of secondary origin particles 
by small bodies. It is in fact an answer to the question 3 
of Zuckerman & Becklin (1993): are the dust grains pri- 
mordial or continually replenished? 

Indeed, wc have shown that if bodies produce dust par- 
ticles with a given size distribution, these particles follow 
very eccentric orbits. Thus, a disk can be seen at large dis- 
tances from the star, and the spatial distribution of dust is 
very close to a power-law. Moreover, the particles have the 
same inclinations as the parent bodies: the thickness of the 
disk thus increases with radius, and the disk is a "wedge" 
disk: (issue 4). The asymmetric spatial distribution can 
be explained if we assume that the parent bodies distri- 
bution is not axisymmetric (issue 2). This assumption is 
not surprising since asymmetry is observed for the FEBs 
and is probably the case in the Solar System for Kuiper 
belt objects trapped in planetary resonances (Jewitt & 
Luu 1995). Since particles are produced with periastrons 
greater than or equal to those of the parent bodies, the 
dust present close to the star is cither produced there or 
brought in by Poynting-Robertson drag: the central part 
must be relatively clear (issue 3). Finally, this model gives 
connection between the particle size and their distances 
from the central star, and enable us to solve the issue of 
the disk mass (issue a). Section 4.4). Thus, the issues 2), 
3), 4) and a) listed in Section 1 are simultaneously solved 
with the simplest form of the model presented here. 

However, in order to explain the flatter gradient ob- 
served in the /3 Pictoris disk than the one obtained with- 
out additional hypothesis, we propose to add two other 
assumptions: 

1) The dust is produced by evaporation of bodies of few 

kilometers in radius, so that there is a balance between 
evaporating gas drag and gravitation of these parent bod- 
ies, and only the smallest particles are extracted. 

2) These bodies are moving slowly towards (3 Pictoris, so 
that they become extinct before arriving close to the star 
where the evaporating rate is large enough to produce the 
largest particles. 

This OEB (Orbiting-Evaporating-Bodies) scenario is pro- 
posed to solve issue 1), but surprisingly it is also able to 
explain the abrupt change in the slope by the cutoff in 
distribution of particle size (issue 5)). It has to be noted 
that the cutoff must not necessary be very sharp and in 
fact is rather smooth in the model since we assumed a par- 
ent body size distribution. The model gives also a natural 
connection between the dust and gas disks (issue b): they 
may well be produced by the same bodies in the same pro- 
cess. The connection between the FEB and asymmetry in 
the disk seems to be possible (issue c)), however this point 
needs further analysis. 



This new model is able for the first time to explain 

simultaneously these issiics, hut wc must now ask if other 
models or hypothesis can also be made with the same 
consequences. If correct, this OEB scenario allows to see 
the P Pictoris disk as a gigantic multi-cometary tail with 
all its components: gas and dust. 
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